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Abstract 

We study two kinds of weighted networks, weighted small- world (WSW) and weighted scale- 
free (WSF). The weight W{j of a link between nodes i and j in the network is defined as the 
product of endpoint node degrees; that is u>ij = (kikj) e . In contrast to adding weights to links 
during networks being constructed, we only consider weights depending on the "popularity" of the 
nodes represented by their connectivity. It was found that the both weighted networks have broad 
distributions on characterization the link weight, vertex strength, and average shortest path length. 
Furthermore, as a survey of the model, the epidemic spreading process in both weighted networks 
was studied based on the standard susceptible-infected (SI) model. The spreading velocity reaches 
a peak very quickly after the infection outbreaks and an exponential decay was found in the long 
time propagation. 
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I. INTRODUCTION 



Complex networks have attracted an increasing interest recently. The main reason is that 
;hey plavan important role in the under standingof complex behaviors in real world networks 



! 



I , including the structure of langua; 
the Internet Q, Q and World Wide Web bl 
chemical reaction networks Il6|], metabolic 



scientific collaboration networks 
12j . power grids food webs LL4 



(it) ], metabolic |l7| and protein networks [l8j . etc. The highly 
heterogeneous topology of these networks is mainly reflected in two characters, the small 
average path lengths among any two nodes (small-world property) Q] and a power law 
distribution (scale-free property), P{k) ~ /c -7 with 2 < 7 < 3, for the probability that any 
node has k connections to other nodes j^j]. 

Most of these networks that have been studied in the physics literature have been binary 
in nature; that is, the edges between vertices are either present or not. Such networks can 
be represented by (0, 1) or binary matrices. A network with N vertices is represented by an 
N x N adjacency matrix W with elements 

1 if i and j are connected, 
otherwise. 



W: 



(1) 



However, many networks are intrinsically weighted, their edges having differing strengths. 
In a social network there may be stronger or weaker social ties between individuals. In a 
metabolic network there may be more or less flux along particular reaction pathways. In 
a food web there may be more or less energy or carbon flow between predator-prey pairs. 
Edge weights in networks have received relatively little attention in the physics literature 
for the excellent reason that in any field one is well advised to look at the simple cases first 
(unweighted networks) before moving on to more complex ones (weighted networks). 

On the other hand, there are many cases where edge weights are known for networks, 
and to ignore them is to throw out a lot of data that, in theory at least, could help us to 
understand these systems better. Recently some weighted network models are proposed. In 



Ref. 



21] , a weighted network model was presented, in which both the structure and the 



connection weights are driven by the connectivity according to the preferential attachment 
rule. In Ref. |22|, a weighted evolving network model with stochastic weight assignments 



was suggested. In Ref. 



it was also found that the strengths of nodes in complex 



networks obey a power-law distribution. In Ref. 



24( , a weighted network was introduced to 



account for possible distinct functional roles played by different nodes and links in complex 
networks by utilizing the concept of betweenness. 

It should be noted that, most of these existing weighted network models are driven entirely 
by the preferential attachment scheme in structure and the preferential strengthening scheme 
in connection weights. In this paper, we will introduce a class of weight to two prototype 
complex networks, the Watts-Strogatz (WS) model and the Barabasi-Albert (BA) model, 
considering link weights depending on the "popularity" of the nodes represented by their 
connectivity. We shall show that the two weighted networks, weighted WS (WWS) and 
weighted BA (WBA), have broad distributions on characterization their structures and new 
results of the spread of infectious diseases. 

II. THE MODEL 

In real systems, a link's weight reflect the familiarity between two nodes, which can 
be regarded as the combination of the nodes' "popularity". For example, in the airport 
network, each given weight is the number of available seats on direct flight connections 
between the airports i and j; in the scientific collaboration network, the nodes are identified 
with authors and the weight depends on the number of coauthored papers. Based on this 
simple mechanism, we will assign weights to links in two prototype complex networks, the 
WS and BA graphs. The two models, which reveal the small-world and scale-free properties 
respectively, have become the most popular models to be used to study the statistics of 
many real networks and explore the dynamics on the network structures. 

WS graph: Starting with a ring of N vertices, each connected to its 2K nearest neighbors 
by undirected edges, and then each local link is visited once with the rewiring probability 
p it is removed and reconnected to a randomly chosen node. Duplicate and self-connected 
edges are forbidden. After the whole sweep of the entire network, a small-world graph is 
constructed with an average connectivity (k) = 2K (in the present work we will consider 
the parameters p = 0.1 and K = 5). 

BA graph: Starting from a small number m of nodes, every time step a new vertex is 
added, with m links that are connected to an old node % with probability 



where kj is the connectivity of the zth node. After iterating this scheme a sufficient number of 
times, we obtain a network composed by iV nodes with connectivity distribution P(k) ~ k~ 3 
and average connectivity (k) = 2m (in the present work we will consider the parameters 
m = m = 3). 

After the networks have been constructed, we shall define the weight w of a link between 
nodes i and j in the network to be the product of endpoint node degrees 

{ (kikj) 9 if i and j are connected, 
= u . (3) 

I otherwise 

where 9 is a positive constant. In contrast to adding weights to links during networks being 
constructed, we only consider weights depending on the nodes connectivity. 



III. STRUCTURAL PROPERTIES 



The standard topological characterization of networks is obtained by the analysis of the 
probability distribution P(k) that a vertex has degree k. Similarly, the first topological 
characterization of weights is obtained by the distribution P(w), which is defined as the 
probability that a randomly selected link between two nodes has a weight value w. In Fig. 
[Tfa)-(b), we plot the cumulative distribution P(x > w/w max ) on linear-log scale for WWS 
networks. The cumulative distribution is defined as 

P(>x) = ]Tp(x'). (4) 

x'=x 

The both distributions decay faster than an exponential decay would. In Fig. ^c)-(d), 
we plot the distribution P{w /w max ) on log- log scale for WBA networks. Fig. [He) shows a 
truncated power law behavior with the apparent exponent a = — 7.10±0.08. Fig. ^d) shows 
a double power law behavior with the apparent exponent a = —2.45 ± 0.03 (for small-w) 
and a = —3.08 ± 0.06 (for large- u>). 

Along with the degree of a node, a very significative measure of the network in terms of 
the actual weights is obtained by looking at the vertex strength s, defined as 

Si = Yl Wi i ( 5 ) 

jeu(i) 

where the sum runs over the set U(i) of neighbors of i. The strength of a node integrates the 
information both about its connectivity and the importance of the weights of its links, and 
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can be considered as the natural generalization of the connectivity. In Fig. Efa)-(b), we plot 
the cumulative distribution P(x > s/ s max ) on linear-log scale for WWS networks. The both 
distributions follow a roughly exponential decay. In Fig. |2^c)-(d), we plot the distribution 
P(w/w max ) on log-log scale for WBA networks. The both distributions follow a power-law 
behavior with the apparent exponent a = —2.66 ± 0.02 (Fig. El^c)) and a = —2.31 ± 0.02 
(Fig. |ld)). 

nn 

The shortest path plays an important role for the transport within a network |25l . I26| . A 
path denotes a sequence of vertices, successive pairs of which are connected via edges. In 
general there exist many paths connecting two given vertices. The shortest path is the one 
with minimum path length among all the paths. In weighted networks, the average shortest 
path length can be defined as 

2 



N(N 



where Wij is the the minimum value of the sum of weights among all the paths between node 
i and j. In Fig. |3] scaling average shortest path lengths of the weighted WS and weighted 
BA networks are shown. The best linear fit in Fig. Efa) yields a slope 0.250 ± 0.001, which 
implies that the behavior of the average shortest path length of the weighted WS network 
remain invariable. Whereas, Fig. Efb) yields that the slope of the linear fit of the weighted 
BA network is about (k) 2e times as that of the SF network. 

IV. EPIDEMIC SPREADING 

In the study of complex networks, a good example is to inspect the effect of their complex 
features on the dynamics of epidemic and disease spreading. It is easy to foresee that 
the characterization and understanding of epidemic dynamics on these networks can find 
immediate applications to a large number of problems, such as computer virus infections, 



distribution of wealth, transmission of public opinion, etc. Recent papers [271 l28l |29[ have 
given some valuable insights of that: for small- world networks, there is a critical threshold 
below which an infection with a spreading rate dies out; on the contrary, for scale- free 
networks, even an infection with a low spreading rate will prevalence the entire population. 
However, so far, studies of epidemic spreading just focus on unweighted networks, and a 
detailed inspection of epidemic spreading in weighted networks is missing. 



To study the dynamics of infectious diseases spreading on weighted networks, we shall 
study the standard susceptible-infected (SI) model j^]. In this model individuals can only 
exist in two different states: susceptible (or healthy) and infected. The model can be 
described in terms of the densities of susceptible and infected individuals, s(t) and i(t), 
respectively, then s(t) + i(t) = 1. Each individual is represented by a vertex of the network 
and the links are the connections between individuals along which the infection may spread. 
In weighted networks, according to Ref. the spreading rate can be defined as 

Xij = ^- (6) 

at which susceptible individual i acquire the infection from the infected neighbor j, where 
w max is the largest value of in the network. In this model, infected individuals remain 
always infective, an approximation that is useful to describe early epidemic stages in which 
no control measures are deployed. 

We start simulations by selecting one vertex randomly and assuming it is infected. The 
disease will spread in the network in according with the rule of Eq. © . In Fig. 0] we plot 
the density of infected individuals versus time in both weighted networks. Since w%3 < 1, 
the smaller 9 is, the more quickly the infection spreads. All the individuals will be infected 
in the limit of long time as lim^oo = 1. We will study the detail of spreading velocity 
at the outbreak moment which is defined as 

V inf (t) = ^*i(t)-i(t-l). (7) 

We account the number of newly infected vertices at each time step and report the spreading 
velocity in Fig. Apparently, the spreading velocity goes up to a peak quickly and leave us 
very short response time to develop control measures. At the moment of infection outbreaks, 
the number of infected individuals is very small, as well as a very long time from the outbreak, 
the number of susceptible individuals is very small. Thus when t is very small (or large), 
the spreading velocity is close to zero. Moreover, in contrast to trivial velocity decay for the 
WWS network (see Fig. Efa)), an exponential decay was found for the WBA network (see 
Fig. 
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V. CONCLUSIONS 



A more complete view of complex networks is provided by the study of the interactions 
defining the links of the system. In real systems, a link's weight reflect the familiarity between 
two nodes, which can be regarded as the combination of the nodes' "popularity" . Based on 
this simple mechanism, we study two kinds of weighted networks, weighted small-world and 
weighted scale-free, in which each link ij in the network has an associated weig 
where ki and kj are the node degrees at the endpoint of the link ij. In contrast to adding 
weights to links during networks being constructed, we only consider weights depending 
on the "popularity" of the nodes represented by their connectivity. It was found that the 
both weighted networks have broad distributions on characterization their link weights and 
vertex strength, an exponential decay for the WWS network and a power law behavior for 
the WBA network. The average shortest path length of the WWS network remain invariable 
comparing with the WS network. Whereas, the average shortest path length of the WBA 
network is about (k) 2e times as that of the BA network. Furthermore, as a survey of the 
model, we investigated the spread of infectious diseases in both weighted networks based on 
the standard SI model. The spreading velocity reaches a peak very quickly after the infection 
outbreaks in both weighted networks. Moreover, in contrast to trivial velocity decay for the 
WWS network, an exponential decay was found for the WBA network. 

As networks play an increasing role in the exploration of complex systems, there is an 
imminent need to understand the interplay between network dynamics and topology. We 
model the weights as solely dependent on the topology, potentially overlooking correlations 
among the weights themselves and the coupled evolution in time. Recently, Barrat et al. pre- 
sented a model .32] to study the dynamical evolution of weights according to the topological 
variations, which is a good starting point towards uncovering the role of such mechanisms. 
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FIG. 1: Linear-log plot of the cumulative distribution P(x > w/w max ) of the scaling individual link 
weight in the WWS network with 8 = 0.5 (a) and 8 = 1 (b) and the log-log plot of the distribution 
P(w/w max ) of the scaling individual link weight in the WBA network with 8 = 0.2 (c) and 8 = 0.5 
(d). The data were averaged on 10 networks of size ./V = 10 6 . 
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FIG. 2: Linear- log plot of the cumulative distribution P(x > s/s max ) of the scaling individual 
link weight in the weighted WS network with 8 = 0.5 (a) and 9 = 1 (b) and the log-log plot of 
the distribution P(s/s max ) of the scaling individual link weight in the weighted BA network with 
9 = 0.2 (c) and 9 = 0.5 (d). The data were averaged on 10 networks of size TV = 10 6 . 
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FIG. 3: Scaling average shortest path length in the WWS (a) and WBA (b) networks. Simulations 
were implemented on the networks averaging over 50 different realizations. 
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FIG. 4: Density of infected individuals versus time in the WWS (a) and WBA (b) networks with 
N = 10 5 . All the plots were averaged over 100 experiments. 
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FIG. 5: Spreading velocity at each time t in the WWS (a) and WBA (b) networks with N = 10 5 . 
All the plots were averaged over 100 experiments. 
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